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Abstract 

This paper shows that the AL (Ablowitz-Ladik) hierarchy of (integrable) equa- 
tions can be explicitly viewed as a hierarchy of commuting flows which: (a) are 
Hamiltonian with respect to both a standard, local Poisson operator J , and a new 
non-local, skew, almost Poisson operator K,, on the appropriate space; (b) can be 
recursively generated from a recursion operator 1Z = ¥LJ~ X . In addition, the proof 
of these facts relies upon two new pivotal resolvent identities which suggest a gen- 
eral method for uncovering bi-Hamiltonian structures for other families of discrete, 
integrable equations. 

Keywords: discrete integrable equations, lattice dynamics, inverse scattering, Poisson geometry, 
bi-Hamiltonian structures. 

1 Introduction 

The Ablowitz-Ladik (AL) system is one of the most studied discrete integrable systems 
of soliton type. It can be thought of as an integrable discretization of the Nonlinear 
Schrodinger (NLS) equation. However, it has recently received a great deal of attention as 
the background lattice system in a variety of modelling applications including optical fiber 
arrays QJE], chaos in dispersive numerical schemes [21H23] . and linkage dynamics |1UII18| . 

AL is an infinite integrable system of soliton type by which one generally means that 
an infinite family of constants of motion for the AL flow can be constructed through the 
inverse scattering transform (1ST) associated to a particular (discrete) eigenvalue problem. 
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In addition, the 1ST framework provides a mechanism through which large families of 
solutions, such as multi-solitons can be explicitly found. 

AL is also an integrable Hamiltonian system with respect to a natural local Poisson 
structure. Indeed this structure is the discretization of the natural Poisson structure for 
NLS. For many integrable partial differential equations (PDEs) of soliton type it has been 
found that these systems are Hamiltonian with respect to two distinct Poisson structures. 
Moreover, these distinct structures are compatible in a way that enables one to construct 
from them a recursion operator which generates the complete hierarchy of commuting 
flows which is the signature of a completely integrable Hamiltonian system and which 
characterizes its Poisson geometry. If a system has two distinct compatible Poisson struc- 
tures it is referred to as bihamiltonian (see Section 2.) In many specific instances the 
bihamiltonian structure provides a means to relate the Poisson geometry to the 1ST via 
resolvent relations associated to the linear eigenvalue problem. 

Examples of bihamiltonian structures related to an 1ST are much rarer for discrete 
integrable systems. The principle goal of this paper is to build and understand this relation 
for AL. 

This paper is organized as follows. In Section 2 we review the second Poisson structure 
for NLS which is due to Magri ^] . Our approach is to explicitly re-derive Magri's structure 
from a Wronskian construction of commuting flows for soliton PDEs due to Calogero and 
Degasperis [7]. We review the derivation of the second structure for the continuous case 
(NLS) because it is an important guide for identifying the second structure in the discrete 
case. A complete and detailed proof that the Magri structure is Poisson may be found in 

In Section 3 the necessary inverse scattering background for the discrete case is pre- 
sented including the derivation, from a generalized Wronskian relation, of the basic op- 
erators £+ and C- in terms of which the recursion operator 1Z and its properties are 
developed. In the last part of this section, the resolvent identities related to these op- 
erators, which are fundamental to the main results of this paper, are presented. These 
resolvent identities represent a novel contribution to the literature on scattering theory for 
discrete systems. Full details of their derivation may be found in |18| . 

Finally in Section 4 the Poisson-geometric interpretation of the recursion operator 1Z 
is made. We prove that AL has a bihamiltonian character (as defined in this section.) To 
establish this it suffices to show that the second structure, /C, for AL is almost Poisson. It 
is an interesting open problem to determine whether or not K. is Poisson. This topic and 
other potential investigations are discussed in the conclusions. 

2 A Generalized Wronskian Approach to the Poisson 
Structure of NLS 

To illustrate the rationale of our approach to the geometry of the AL equations ^Hj , we 
first look at the parallel continuous object (the NLS hierarchy), where the geometry is 
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already understood. 

The NLS equation Q is a well-known integrable PDE. As such, it possess, infinite 
families of linearly independent constants of motion in involution, and families of explicit 
special solutions such as N-soliton solutions. In the late 70's, Magri showed that NLS has 
a bi-Hamiltonian nature |19| . This property (which now characterizes many integrable 
PDEs), means that the NLS equation can be written as a Hamiltonian system with re- 
spect to two different, independent Poisson brackets. By composing the Poisson operators 
(denoted J and K, for instance) induced by these brackets in the appropriate manner, one 
obtains a recursion operator, R, capable of generating a commuting family of Hamiltonian 
flows, which include NLS. Typically one of the aforementioned Poisson operators, say J, 
is invertible and then R = KJ~ X . 

The central idea is to recover Magri's Poisson structure for NLS from two integro- 
diffcrcntial operators associated to the AKNS hierarchy, a collection of commuting inte- 
grable evolution equations which includes NLS. 

The first of these operators, which we denote L+, was constructed by Calogero and 
Degasperis using a generalized Wronskian technique |7|. The basic idea behind the gen- 
eralized Wronskian approach to integrable evolution equations, is to generate a set of 
scattering relations between the asymptotic behaviors (as x — > +00 and x — > -co) of the 
wave-function solutions of some initial eigenvalue problem, such as the Zakharov-Shabat 
problem 



The second operator, L_, appearing in Calogero-Degasperis' work and giving rise to 
the AKNS hierarchy is related to L + in a sort of "adjoint" way as described below. 

The specific form of the operators L + and L_ (which act on a certain space C of rapidly 
decaying complex-valued functions (q(x), r(x)) T -the potentials,) suggests the construction 
of an anti-symmetric operator K arising from L + , L_. 

As shown in ^B] (see below) one can indeed use L+ and L- to define two geometrically 
meaningful operators K and R which act on C . The first one is Poisson, and gives rise to 
Magri's bracket for NLS. The second one is a recursion operator for the AKNS hierarchy 
which may in fact be portrayed as the composition of K with the standard Poison structure 
for AKNS (given essentially as multiplication by i). 

The specific structure and relationship between L + and L_ in the continuous setting 
will eventually guide the construction of their discrete counterparts C + and £_ and also 
motivate the definition of K, and 7Z -discrete analogs of K and R for the AL hierarchy 
[T8| . These results are presented next. 

2.1 Obtaining Magri's Poisson structure, K, for NLS through L + 



In their 1976 paper, Calogero and Degasperis begin with the the eigenvalue problem l(T|l 
and use a generalized Wronskian technique to arrive at the following formally defined class 




(1) 



and L 
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of integrable equations: 



(r t {x, t), -q t (x, t)f = l(L+)(r(x, t),q(x, t)) 1 



(2) 



where 7 is an entire function of the integro-differential operator 



+ 2i 



rl+q 



-tIa-T 



ql+q -ql+r 



U = 



r+oa 

/ ( • )dy. 



(3) 



Perhaps one of the best-known procedures for generating hierarchies of non-linear in- 
tegrable evolution equations is due to the work of Ablowitz, Kaup, Newell and Segur 
(AKNS) 1 . In our context of interest, their method utilizes once again the Zakharov- 
Shabat eigenvalue problem together with a second linear operator prescribing the time 
evolution of the wavefunctions. The hierarchy of integrable evolution equations arises then 
as a series of compatibility conditions associated to the linear problems just described. All 
the evolution equations in the hierarchy stem from the same eigenvalue problem, yet each 
corresponds to a different time-evolution operator. 

Calogero and Degasperis write the AKNS hierarchy as 



(r t (x,t),-q t (x,t)) T = -2A{LJ){r{x,t),q{x,t))' x 



where 



1 

2i 



1 

-1 



d x 



rl-q 
ql-q 



-rl- 

-ql- 



7_ 



)dy. 



(A is an arbitrary entire function), and show that J5J and Q are equivalent. 
The proof of this fact is entails showing that 

Ll((r : q) T )=L n _((r,qf), 



(4) 



(5) 



(6) 



which follows, in turn, from the observation that L!?((r, q) T ) is in the kernel of the differ- 
ence operator L_ — L + -the continuous "kernel condition." This can be verified directly 
for n = and other low- values of n. The proof for general n (which relies on the definition 
of two auxiliary functions of the complex variable z) can be found in the Appendix of the 
original Calogero-Degasperis paper. 

We remark that, under the identification r = —q, the system of coupled integrable 
evolution equations reduces to a system of evolution equations for a single field q, 
comprising the NLS hierarchy, which contains the well-known NLS equation 



-iqt = q xx + 2|<?|V 



(7) 



In other words, in the reduction r = —q, (0J becomes the NLS family of equations. 

Let us now consider the space of complex, vector valued functions of a real variable x 
given by 

C = {(q(x),r(x)) T 6 C 2 : \q\, \r\ -> as x -> ±00} . 
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The tangent bundle to C may be endowed with a non-degenerate bilinear form locally 
described by the inner product 



((vi,v 2 ) T ,(w 1 ,w 2 ) T ) = - j v 1 w 2 + v 2 wi dx, (8) 

where, (v\, v 2 ) T , (u>i, w 2 ) T are tangent vectors to C at the point (q,r) T , and J f dx = 
/ dx. One may also define a skew-adjoint operator J = Diag(— on the tan- 
gent bundle to C which, together with the inner product just defined, gives rise to the 
complexified standard Poisson bracket on C , that is, 

{F,G}j = (VF,Diag(-M)VG) = ({-S r F, ~6 q F) T , DiagH, i)(-S r G, -5 q G) T ) . 

Here, F and G are functionals on C which become real-valued in the reduction rk = qk 
defining the standard setting for NLS. This reality condition implies that F and G must 
be symmetric (or anti-symmetric) in r and q. Furthermore, (5 q ( ■ ),5 r ( ■ )) T denotes the 
variational derivative, whereas V( ■ ) = (— S r ( ■ ), — 6 q ( ■ )) T gives the functional gradient 
of each such functional, with respect to the given inner product. 

We now re-consider the operators L + and L_ in the context of the function space C. 

The presence of the integral operators I + and /_ suggests that L + and L_ may be 
combined so as to yield an anti- symmetric operator, K related to the hierarchy. Indeed, 
as we will show below, the special properties of L_ + L + and L_ — L + allows us to recover 
both: the (second) Poisson structure for NLS discovered by Magri 19;. (prescribed by the 
Poisson operator K), and the recursion operator R = KJ^ 1 associated to the hierarchy 
of flows, in the current setting. We begin by defining, 

\0 1) \r(I--I+)r -r{I--I+)qJ w 

This operator K, just as J, acts on the tangent space to C at a point and defines a 
(point-dependent) operator on the tangent bundle of C . 

Proposition 1 The bracket {F, G}k = (Vi 71 , K^JG) associated to K is skew- symmetric 
and satisfies the Jacobi identity. Hence, K is Poisson. □ 

A proof by direct-calculation of Proposition 2] can be found in ^SJ. The technique relies 
on a standard integration by parts formula (see Appendix A in [T5J , for instance.) 

There are at least two possible avenues to establish a precise link between the Poisson 
operator K and the hierarchy of integrable evolution equations (QJ. The first one begins 
with the observation that, in the reduction r = —q, Hamilton's equations associated to 
K become those given by Magri's Poisson bracket ^Hj- Since the integrable hierarchy in 
the previous section is the NLS hierarchy for r = — q, and Magri's bracket is known to 
prescribe a bi-Hamiltonian structure for NLS (together with the Poisson bracket given as 
multiplication by —i), the above observation connects K with the AKNS hierarchy via the 
identification r = —q. 
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A second and perhaps more constructive approach relies on directly checking that K 
and J form a compatible pair of Poisson structures and that the recursion operator KJ^ 1 
does indeed generate the AKNS hierarchy of flows. This approach has the advantage of 
giving an explicit description of the bi-hamiltonian nature of the evolution equations in 
the AKNS hierarchy. 

Recall that the compatibility of two Poisson operators such as K and J amounts to 
proving that their sum is also a Poisson operator on the manifold in question. Two 
compatible Poisson operators give rise to a bi-Hamiltonian system when there exists a 
vector field that is Hamiltonian with respect to both Poisson structures I n this case, 
one can generate a hierarchy of bi-Hamiltonian vector fields, by recursively applying the 
composition of one of the Poisson operators with the inverse of the other to the original 
bi-Hamiltonian field (assuming, of course, that one of the Poisson operators is indeed 
invertible.) We define 



R = Kr 



1 
1 

1 

-1 



d x 



-q(I_ - I + )r q(I_ - I + )q 
r(I- - I+)r -r(/_ - I+)q 



i 
-i 



(L_+L, 



-1 

1 



(10) 



and consider the functional H = — J rq. Then, J\7H = Diag(— i, i)(q, r) T = (—iq,ir) T , 
and X n = R n J\7H defines a hierarchy of bi-Hamiltonian vector fields. Hamilton's equa- 
tions associated to these fields are 



((q,r) T ) t = R n JVH. 



(11) 



Notice that X = JVH = (-iq,ir) T and Xi = KVH = {q x ,r x ) T . The standard NLS 
flow is prescribed by 



^3 = i [(&*, -r xx f + (-2rq 2 ,2qr 2 ) 



2\T1 



(12) 



The hierarchy {TTJ is the NLS Hierarchy. 

We shall next show that the NLS hierarchy is equivalent to a sub-hierarchy of 
the family of Calogero-Degasperis flows In fact, we will demonstrate that the entire 
Calogero-Degasperis hierarchy can be generated from polynomials in R. The argument's 
basic idea is to write the recursion operator R in terms of L_ + L + and exploit the 
relationship between L + and L_ given in ©. 

Proposition 2 The NLS Hierarchy, obtained by recursively applying the recursion op- 
erator R to the (Hamiltonian) field (~iq,ir) T , is equivalent to a sub-hierarchy of the 
Calogero-Degasperis family J2J. 

Proof. We re- write Hamilton's equations 1)11(1 . using the form of the recursion operator 
given in ffTTTfl . to obtain 



((q,r) T ) t = a 2 (L_+L + T(r,q)' 



(13) 



G 



where 



02 = 



c 



-i 

1 



) 



(14) 




(15) 




(16) 



(17) 



which is clearly a sub-hierarchy of the Calogero-Degasperis family of flows. 



□ 



Using Weierstrass approximation one can replace the monomial in 117fl by an arbitrary 
analytic function, A, of L + and finally re- write the equation as ((r, — q) T ) t +A(L + )(r, q) T = 
0, which is the form of the flow derived by Calogero and Degasperis. Hence, 

Corollary 1 The entire Calogero-Degasperis family of flows may be formally generated 



3 Inverse Scattering Preliminaries 

This section begins with the direct scattering formulation for the Zakharov-Shabat eigen- 
value problem (|18fl originally studied by Ablowitz and Ladik in connection with the inte- 
grable equations bearing their name. 

After identifying (and characterizing) two pairs of eigenfunctions (f>k(z), <f>k{z), tpk(z), 
tpk{z) (which behave nicely at either — oo or oo), we focus on a Wronskian relation intro- 
duced by Ladik and Chiu as a means to study a family of integrable evolution equations 
arising from the linear problem i|18[) (which includes AL.) 

The generalized Wronskian identity presented by Ladik and Chiu provides, once again, 
a way to relate the evolution of the potentials and (defining (|18fl ) , and the kinematics 
of the scattering coefficients (i.e., the parameters specifying the asymptotic behavior of 
the eigenfunctions mentioned above.) By so doing, it also singles out two sum- difference 
operators L and L~ l , which together generate the collection of (integrable) equations 
prescribing the potentials' evolution. 

As we will soon see, the specific form and relationship between L and L~ 1 will guide 
the construction of two new operators £+ and associated to a sub-hierarchy of the 
Chiu-Ladik equations which encompasses the AL system |18j . The latter operators can 
be viewed as discrete analogs of L + and L_ in more than one way. In particular, they 



from polynomials in R. 



□ 
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will play a key role in unveiling the Poisson-geometric picture behind AL, as we begin to 
indicate below. 

Exploiting the time-independence of one of the scattering coefficients (an essential 
feature of these kind of integrable equations), one can obtain both: a sequence of constants 
of motion, and a generating function for the associated hierarchy of variation (gradient) 
fields. It turns out that it is also possible to explicitly write down (resolvent-type) identities 
giving the infinite hierarchy of gradient fields in terms of powers of L and L~ x . As explained 
below, these identities will be instrumental for exhibiting the bi-Hamiltonian nature of the 
AL hierarchy in Section 0] 



3.1 A generalized Wronskian leading to AL flows 

In 1976, Chiu and Ladik orchestrated a generalized Wronskian identity based on the 
discrete version of the Zakharov-Shabat eigenvalue problem 

which generated a large family of integrable evolution equations associated to AL. Here, 
z £ C is the eigenvalue parameter, and the complex-valued potentials q k and r k are 
assumed to vanish rapidly for k — > ±oo. 

Let $k(z) and &u{z) respectively denote matrix solutions to (|18fl and to a second eigen- 
value problem of the same form but for a new set of potentials q' k , r' k . The initial Wronskian 
identity proposed by Chiu and Ladik relates <& k and <&' k to a,a,b,b and a',a',b',b (the 
"unprimed" and "primed" scattering parameters) via the right-hand-side of the equation 
below 0: 



P(k + 1) M(k,z) $ fe+ i = & k T P(k)F(k)$>k\°?oo' 

k—~~oo 

l/z -q k 



where M(k, z) = 



q< k l/ Z J F ^ k + 1 ^ F ^\-r k 



(19) 



Above, P(k) = YiJLki^ ~ r kQ_k) and F(k, z) is an arbitrarily chosen matrix. The validity of 
(|19l) results from the telescoping of the terms in the specified series, which can be checked 
directly (see ^B] for details.) 

Due to the asymptotic decay of the potentials, the right-hand-side of Ijl9(l makes sense 
and may be written in terms of the scattering parameters (for the "primed" and "un- 
primed" linear problems.) Its precise form depends both on the choice of matrix F{k, z) 
and on the particular eigenfunctions making up the columns of & k (z). 

Remark 1 We will use the notation f k and f(k) interchangeably, depending on emphasis, 
to denote functions of the discrete variable k. The z dependence of objects such as <fr k , 
F(k) and the scattering parameters a, a, b, b, may or may not be explicitly noted, depending 
on the context. 



Following the work of Chiu-Ladik, we now specify the precise form of the right-hand-side 
of (|19|l by choosing a particular fundamental matrix $> k (z) and determining its asymptotic 
form. We begin by considering the four special solutions 



<Pk{z) ~ z fc (l,0) T , $ h (z) ~ z~ k (0, 1) T , for k -» -oo 
Mz) ~ z~ k (0, if, Mz) ~ z k (l, Of, for k -> oo. 
of 118(1 and the scattering relations prescribed by 



4>k(z) 
<i>k(z) 



b(z) a(z)\ (Mz) 
a(z) b(z)J {Mz) 



= S(z) 



We then set 3>fc(z) = (<fi k (z), ^fc(^)) and determine that 



o ^V fc 



for fc 



4>a 



>z~ 



ipk(z) 







62 ' 



-k ! fOT fe 



(20) 
(21) 

(22) 
(23) 



where det£(z) = aa — bb = Co, and Co = lim / r.^_ oc P(fc) = irff^l — r k q k ), as shown by 
Ablowitz, Prinari and Trubatch 0. (See also ^Hj for more details.) 

We now adapt the results established by Chiu and Ladik to our choice of &(z) and 
arrive at two Wronskian-type identities based on lfT§jl . They are: 



E P (k + lWk 



k= — oo 



2A l 



^ 0) (fc) N 



fe+i 



f alb a! — a\ 
ab' 







C' 



E z 



(*)**)!-»> (24 ) 



E p ( fc + 



k= — 00 



J(A^) 



-W(hf\k) 



fe+i 



1 

37 



r/i/ 



OO 
fcj I— OO 5 



(25) 



where the integro-differential operators A, A 1 are given by 



A 







-r k J k { q] E-) r k J k {r' J E+) 
-q' k Jk(q 3 E-) q' k J k (r' j E+) 





q' k q k E~ 









-E+iqkJkiq'j)) E+(q k J k (rj)) 



(26) 



9 



*-i = (E + \ (r' k J k {q' J E+) - r ' k J k ( rj E-)\ (-r' k q> k E+ r' k r k E- 
E-J + \q k JkWjE + ) -qMriE-))^\ 



+ (1 - r k q k ) 





E-(q' k )Jk(Qj) -E-(q' k )J k (r' 3 ) 



(l-rU k )( E+ir ^ M) ~ E+ ^m, (27) 



as demonstrated in the original Chiu-Ladik work 8 . The shift operators E^ act in the 
standard fashion, E ± v k = v k ±i, whereas J k (v(j)) u = J2jLk v U) u (j)' and Jk(v(j)) u = 
p(k)/2 + J2f =k [(Uizl(l - r&)/(l - r m ))(v(j)u(j))/(l - r jqj )}, for p(k) = pll£ fc (l ~ 
r fc9fc)/(l — r kQk) and p an arbitrary constant. 

As explained both in ^8] and in Chiu-Ladik's paper, Eq. I|24(l arises from an iteration 
in powers of z 2 prescribed essentially through a suitably chosen ansatz for the M(k, z) 
factor in the summand of 11911 . Specifically, 



_ /y 2fr(l-i)^M _i_ ft(0 
where 

We remark that such an ansatz completely determines the structure of z). The second 
equation emerges in a similar manner from an iteration in powers of 1/z 2 . 

The pivotal identity for obtaining the precise time-evolution of the vector potential 
{r k ,q k ) and the scattering parameters associated to (I18fl . arises by taking linear combi- 
nations of (|24[1 and 11' ">D for different values of I > 18 . When temporal considerations 
are brought to the forefront, the aforementioned combination of (|24() and (|25() becomes 
time-dependent (through the ^-dependent vector potential, scattering data, etc.) One 
then exploits this time dependence by defining the scattering parameters in the "primed" 
eigenvalue problem as time-evolved variants of those in the "non-primed" scenario, namely, 

a' = a(t), a' = a(t), b' = b(t), V = b(t), 

r' k = r k (t), q' k = q k (t), $' fe = C = C (t). (30) 

As detailed in Chapter 3 of JSj , evaluating the right difference-quotients and limits -based 
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on the substitutions i|3U[l - leads to the operators 
(E~ 0\ f-pE-(r k )\ 



+ (1 - nt<7fc) 



+ (1 - r k q k 




-r k J+(q 3 E-) r k J+( rj E+)\^f 



+ \~q k J+( qj E-) q k J+(r]E+)) + {qlE- -q k r k E+ ) ' (31) 



o 2H +(1 - r * ft) U-( ft ) 








V 1-rjQi J 1 






\. 1 - r j1j)j 







r^) E~ 




r k J+( qj E+) - rk J+( rj E-)\ f-r k q k E+ y\E~ 
q k j£{q E+) -q k Jt{r 3 E-))^\ 



(32) 



(arising from A and A 1 respectively), and to the following (simple!) time-evolution 
equations for the scattering parameters: 

((a(z),a(z)f) t = 0, ({-b{z)Mz)) T ) t =uj(z 2 )(b(z)Mz)) T , (C )t = 0. (33) 



Remark 2 The operator present in (|3 II) and 132fl is defined by j£ (v,j ) = Y2°°- k u j 
and is the discrete analog of the integral operator i+ , appearing in the expression for the 
continuous integro- differential operator L+ given in 0. The symbols p^p denote (discrete) 
integration constants. 

Obtaining evolution equations for the potentials r k {t), q k (t), requires more careful argu- 
ments which exploit the analytic properties of the Jost functions associated to the linear 
problem l|18|l and its four special solutions (|2U[1 - H21[1 . Through the proof of Theorem 1 in 
|18|. one arrives at the following time-flows: 

((-r k ,q k ) T ) t =uj(L)(r k ,q k ) T , (34) 
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where uj(x) = uii(x) + w 2 {x~ 1 ) and L and L^ 1 are the sum- difference operators above. 

The family of discrete integrable evolution equations (|34[l arising from the Chiu-Ladik 
Wronskian approach includes several discrete versions of mKdV equations as well as the 
well-known discretization of NLS constructed by Ablowitz and Ladik, namely, 

i{Qk)t = (1 + |<?fc| 2 )(<7fc+i + Qk-i), (35) 

(up to the linear term —2q k .) The latter -which we will also call the AL equation- can 
essentially be obtained from one of the simplest symmetric w(£), namely L + L^ 1 . Indeed, 
as one can directly check, 

((-r fe , q k ) T ) t = -i(L + L- 1 ){r k ,q k ) T = - r k q k )(E+ + E-)(r k ,q k ) T (36) 

In the reduction r k = —q k , Eq. becomes . 

3.2 The constant of motion log a and associated gradient fields 

Recall that the scattering parameter a is an analytic function of the complex variable 
z on the complement of the unit disc, and can therefore be represented by its Laurent 
series expansion (see |18| and [3], for more details.) Since a is time-independent -as made 
evident through the evolution equations l|33(l - so are the coefficients of its Laurent series. 
Such coefficients then constitute an infinite family of constants of motion associated with 
the flows B34JI . The same reasoning shows that the coefficients of the Taylor series for a 
(on \z\ < 1) also give an infinite family of integrals. 

In their 1975 paper, Ablowitz and Ladik utilize a recursive technique to arrive at 
the collection of constants of motion arising from log a: 

Ci=^2q k r k+1 , C 2 = ^V fe+ ig fc -i(l - r k q k ) - ^q 2 k r 2 k+1 , ■■■ (37) 

The same formal method (outlined in |18p. may be used to compute the hierarchy of 
integrals associated to log a, namely, 

Ci = ^r k q k+x , C 2 = ^2 q k +ir k -i(l-r k q k ) - -rlql +1 , ■■■ (38) 

One may also obtain a hierarchy of "symmetric" constants of motion, by combining the 
above as C , C x +C u C 2 +C 2 , where C = UZ-oci 1 - r kq k )- (See £3, E3-) 

We now determine a hierarchy of gradient fields associated to the conserved quantities 
arising from the series expansion of log a and log a just discussed. Our approach here draws 
on the standard treatment of the parallel derivation for the modified Korteweg-DeVries 
equation (mKdV), a well-known integrable PDE on the real line |20| . 

The argument is based upon considering a variation (q k ,f k ) T = (a k ,b k ) T of (q,r) T 
along the space of rapidly vanishing (vector) potentials, while keeping the eigenvalue z 
fixed. By finding the expression for the induced variation on the special eigenfunctions 
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<f>k:^>k,ipk:^k (the special solutions given in (|2U[I. (|2T|) to the linear problem (EHI)j one 
is able to precisely calculate the effect of the potentials' perturbation on the scattering 
parameters a, a, b, b. For instance, 



leads to 



^-^nf^fc + ir 1 ^ «j) ^ (39) 



(40) 

in the manner outlined in Appendix 14.31 of this work. This, in turn, allows us to deduce 
an explicit expression for the variational derivative o^-loga = (5 log a/Sq k , S log a/Sr k ) of 
loga, namely, 

S k (loga - £ log(l - Pfcft )) = ^^(^V«^T • (41) 

By Laurent-expanding the expression for this generating function one can then obtain 
the desired hierarchy of "gradient" fields, each given by the variational derivative of the 
corresponding coefficient in the expansion of loga. 

The results just outlined, together with the supporting arguments presented in Ap- 
pendix prove the following theorem. 

Theorem 1 The (z-dependent) gradients o/loga(z) andloga(z) are generating functions 
for the hierarchy of variation fields 8^C n and S k C n associated to the constants of motion 
(|38|l and l|37l) respectively. Their explicit expressions are: 

6 k lo g a + S k H = ^ll(4 2 i^ 2) , -t&tPr, (42) 
S k \oga + 6 k H = ^±^(-4>^\ ^IM ] ) T , (43) 

where H = -YlkL-oo^gi 1 ~ r k1k) = -logC and P(k + \) = ]\ c * =k+1 (\ - r k q k ). □ 

Notice that, on the unit circle \z\ = 1, both Eq. i|42|) and Eq. (|43|l make sense (as the 
scattering coefficients a{z),a{z) are holomorphic on the regions \z\ > 1 and \z\ < 1 of C 
respectively.) One could then combine them to obtain a "symmetric" representation of the 
variational derivatives. This line of reasoning leads to the symmetric constants of motion 
stemming from the combination of (|38|l and (|37|l . 

3.3 Constructing the generating operators: £ + , £ and 1Z 

Prima facie, one observes several qualitative analogies between the operators (defined 
in Eq. associated to continuous NLS and L + L -1 . Both are second order operators, 
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and both may be obtained via parallel (generalized Wronskian) techniques (in the discrete 
and continuous settings respectively.) Furthermore, they each give rise to (continuous and 
discrete) integrable versions of NLS in a parallel manner: 

i((-r, q) T ) t = L 2 + (r, q) T = -1/2 {{l/2)r xx - r 2 q, (l/2)q xx - q 2 r) T , 
i{(-r k ,q k ) T ) t = (L + L~ 1 )(r k ,qk) T = (1 - r k q k ){E+ + E~)(r k , q k f . 

(Note that the later equation is the differential-difference version of NLS (up to a linear 
term) arising in the work of Chiu and Ladik, just presented in Eq. (j5"B*|> in the previous 
section. As pointed out earlier, we will also refer to this equation as the AL equation.) 
More precisely, setting £ + = L + L^ 1 , we note that, based on Eqs. (|31fl . (|32|) . 



C+ = (E+ + E- 



' + (1 
+ (1 - r k q k 

+ (!• 

+ ( r kJ k (qj) 



»•*?*) [-p {E 




q\E -q k r k E^ 



(44) 



Remark 3 From this point on, we will assume that the integration constants p, p vanish, 
unless explicitly noted otherwise. 

Proposition 3 The operator C + generates a sub-hierarchy of the family of evolution equa- 
tions l|34|) given by 

((-r kl q k ) T ) t =V(£+)(r k ,q k ) T , (45) 
where V denotes an arbitrary polynomial in C+ . 

Proof. Since L and commute, all powers of = L + LT 1 can also be written as 
linear combinations of powers of L and L -1 , and hence P(£+) = lo\(L) + ^(i -1 )- □ 

Note that the hierarchy (|45|l encompasses various versions of the AL equations. In the 
reduction r = —q, these include (|35|l . and also the more standard 

i(q k ) t = (1 + \q k \ 2 ){q k +i + q k -i) - 2q k , (46) 

obtained from P(C+) = £+ — 2. Observe also the parallel with the Calogero-Degasperis 
hierarchy (J2J associated with NLS. 
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In the continuous setting, we explored a second recursive-type operator giving rise 
to the NLS hierarchy generated by L + . This operator, denoted by R, was essentially 
constructed as the sum of L + and L_ (see (0 and (0) and was later seen to carry 
important geometric information concerning the Poisson geometry of NLS. 

Our next goal is to construct discrete analogs of the operators L_ and _R, and to exhibit 
their connection to the AL hierarchy (|45fl stemming from C + . The precise geometric 
implications of these constructions will be addressed in Section 

We begin by observing that the form of L_ (specifically the presence of the integral 
operator i_ in place of the /+ appearing in L + ) suggests an adjoint-type relation between 
L + and L_. 

Proposition 4 Let 03 = Diag(l, — 1), and let * denote the adjoint of a given operator 
with respect to the inner product ©. Then i_ = a^L^a^ 1 . 

Proof. We first compute L* + using the aforementioned inner product and integration by 
parts. The proposition then follows upon conjugating the calculation's result by 173 |18j . 
□ 

In the remaining part of this section we construct the adjoint of C+ and use it to 
define the discrete counterpart of L_. We then show that, once again, the sum 
C + + C- (appropriately conjugated) defines a (discrete) recursion operator 1Z for AL 
which is analogous to the "continuous" R associated to NLS. 

Let (■, ■) denote the (non-Hermitean) inner product on the complex space • ■ ■ C 2 x C 2 x 
C 2 • • • (indexed by k in Z) defined by 

((a, 6) T , (c, df) = -J2 [a k d k + b k c k ) (1 - q k r k )'\ (47) 

where (a, b) T = (a k ,b k ) T , (c, d) T = {c k ,d k ) T , k e Z. 
Inspired by the preceding proposition we define 

= a 3 C* + a^ = a 3 (L* + (Z^ 1 )*)^ 1 , (48) 

where * denotes the adjoint of an operator acting on ■ • • C 2 x C 2 x C 2 • • ■ with respect to 
the inner product l|47|l just defined. We then have, 

L* = (1 - r k q k ) 
+ (1 - r k q k 

+ (1 - r k q k 




r k q k 

r k J k (q,E-) r k J-{r E+)\ [0 -r\E+ 
-q k J k (q 3 E~) -q k J-( rj E+)J ^ \0 q k r k E+ 



(49) 
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v 1—rkQk > 



-r k J~( qj E+) -r k J-{r 3 E-)\ (r k q k E+ 
q k J k ( qj E+) q k J-{r 3 E-) )^\-qlE+ 



(50) 



where J k denotes the equivalent of the integral operator 7_ in the discrete context, that 



is, J k (uj) = X^ 7 -=-oo u j- Summing the two expressions above and conjugating by 03, yields: 

— ' ~ r k J k \T=^) 

1xJ k (1^7) 



->3 

C- = (E + + E~) - (1 - r k q k ){E+ - E~) ' " ' v ' ' '"' > " v ' - •</• 



(1 - r fcgfe )£+ I 

\ \-r k q k 




(r k J k {qj) r k J k ( rj )\ . + _ , _ /r fc g fc ^ -r 2 k E+\ 



The (lengthy but straight-forward) calculations yielding formulas (|49(l for L* and (|50|) for 
(L^ 1 )* are outlined in Appendix B of |18| . 
Let us now introduce the operator 



K = 




( H ^ 






') 


r:r< 









+ (1 - r k q k )(E+ + E~) \ ^ r r f k 

• l-r k q k l-r k q k , 

Ik (4 J k ) fa ) -q k {4 J k ) fe )\ 
~n (4 - J k ) {r 3 ) r k (J+ - J' ) {q 3 ) ) 



-Qkn -qi 

-rl -r k q k 



(E++E-). (52) 



1G 



Based solely on its form and method of construction, 1Z may be regarded as a discrete 
analog of R. The analogy is strengthened by the fact that recursive applications of TZ 
to the (Hamiltonian) field (iqk, —ifk) T yield the same hierarchy of evolution equations 
generated by £ + . This can be verified by direct calculation for the first few iterations; it 
will be proven in the general case in Theorem 0] For the moment, note -once again- the 
parallel with the continuous setting. 

In Section 0] we will also prove that TZ is in fact the true geometric analog of the 
recursion operator R associated to the the NLS hierarchy. In particular, we will realize TZ 
as the composition of two skew-symmetric operators J~ l and IC on the appropriate space 
of complex, vector valued functions of a discrete real variable fc. 

The next section contains a pivotal result for rigorously proving the geometric character 
of 1Z; thus the framework for unveiling the Poisson geometry of the AL equations begins 
to emerge. 

3.4 Resolvent identities for L and L 1 

In this section we establish two resolvent-type identities. The first one links the hierarchy 
of evolution equations generated by powers of L to the hierarchy of gradient fields (|42|) 
stemming from the generating function 8k log a + SkHg. The second one, relates the flows 
given in terms of powers of L^ 1 to the hierarchy of gradient fields l|43[) stemming from the 
generating function 5k log a + 8k Hq. 

The fact that the operators L and L _1 satisfy these identities will be key for postulating 
and proving several fundamental results in Section 0] . The first of these is Theorem 
followed by Corollary [3| and Corollary Through them, we show: (1) that the hierarchy 
of flows (|65|l generated by the recursion operator 1Z is equivalent to the hierarchy arising 
from powers of C+ = L + L^ 1 ; and (2) that the Calogero-Degasperis scheme leading to Eq. 
© (the continuous "kernel condition") and Proposition^ extends to the discrete setting. 
Furthermore, once Theorem 0] is established, we use the resolvent formulas to prove that 
the ^-generated flows are bi-Hamiltonian (Theorem [5]) 

In the study of hierarchies of integrable evolution equations, it is morally expected that 
the operator playing the role of L or L _1 (or even L + in the continuous case,) will satisfy 
an identity with the properties just described. This is because the generating function 
(Jfcloget (or its equivalent) comprises a full hierarchy of gradient fields which we know 
should be recursively obtainable from L (or its corresponding counterpart.) In essence, 
the resolvent identity allows us to compactly write the same hierarchy of gradient fields 
but now in terms of L and its powers. 

Based on this line of reasoning, one expects that the resolvent acting on the simplest, 
initial flow in the L-hierarchy (or L~ ^hierarchy) should yield the full sequence of gradients. 
It turns out that, indeed, 

Theorem 2 Let L, L^ 1 be as in (|32fl respectively. Then L and L 1 satisfy the 
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following resolvent identities 



(I - z 2 L) 1 (r k ,q k ) T = (1 - r k q k )(8 k \oga + 5 k H ), 
(I - z 2 L- 1 y 1 {r k ,q k ) T = (1 - r k q k ){S k logo + 5 k H ), 



(53) 
(54) 



where H = - J2T=-oo lo g(! _ 



Proof. Setting (/3 k , -a k ) T = P{k + l)/a (V'fVfc+i, -0fc 1) V'fc+i) T , in accordance with © 
in Appendix II, we recall from Theorem that the non-standard squared eigenfunctions 
given by —a k and /3 k encode the hierarchy of gradient fields stemming from the generating 
function S k loga+5 k Ho. The resolvent formula i|53|) we seek to establish explicitly portrays 
(— a k ,/3 k ) T as a generating function for the hierarchy of fields arising from powers of L. 

Guided by the form of the leading term of L we select the first component of (|101|) and 
the second component of (|100|) (re- writing the 07-factor in each using Hl(J4|i ) and combine 
them so as to get 



E- 






E+ 



(1 - r k q k ) 



ftk 



-a k 



(1 - r k q k ) 



Pk 
-a k 



-(r fe _i +z 2 r k ) 







E- 



J k+i(lj) - J k+i( r j) 



(3k 
-a k 



q k+1 + z 2 q k J \ 

Re-grouping terms and dividing through by — z 2 , the last equation is equivalent to 

Pk 



I-z' 



E~ 






E^ 



(1 - r k q k ) 



~OL k 



+ z- z B = 



(55) 



(56) 



where B is the second term on the right-hand-side of (|55|l . 

Based on the expected form of the resolvent identity and the formula for L, Eq. (|56() 
suggests that B ought to simply amount to all but the first term of the operator — L acting 
on the vector [(1 — r k q k )(/3 k , —a k ) T \. In other words, 



o 

E^ 



+ 



-Tk-l 









(1 - r k q k ) 
' Pk 



(1 - r k q k )T{k) 
Qk+lJ \-a k 



Pk 
T(fc) 



- J k(lj) J k( r j) \ ,r.-, 

4+M -■tf+ifa), 1. 



(where the symbol = signifies a claim yet-to-be-established.) Expressing B (in (|55fl) in 
terms of T(fc), one notices that (|57|) is equivalent to Diag(— r k , q k )A = 0, where 



A = 



' rk - iqk n +z2 2 )r(k)-?(k)( E ~ 



+ (1 - r k q k ) 



Pk X 
-a k . 



, (58) 
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which in turn holds if A = 0. Iterating once, we re-write the factor Diag(i?~, E + )[(l — 
r kqk)(f3k, —o>k) T ] in the second term of A with its equivalent form given in Eq. J55j). 
Grouping terms in z 2 , we see that A amounts to the following quadratic equation in T(fc): 



rk-iqk 
qk+irk 



T(fc) 



z 2 T{k) 











Qj 



T(fc) 





Qj+i 
( QjPj+rjaj 



Pk 



(59) 



Now, due to the form of the matrix T(fc), a simple calculation shows that the coefficient of 
the z 2 term in Eq. I|59|l vanishes. The quadratic term within the first line of 159fl telescopes 
down to its boundary term, J k (qj(3j + rjOtj)(qkrk-i, —qk+irk) T , which is evidently the 
negative of the linear term in the first line of (|59|l . Hence A = and the resolvent formula 
follows. □ 



4 Geometry: almost Poisson structure of AL 

In this section we exhibit the bi-Hamiltonian character of the AL hierarchy of equations, 
as defined in Eq. 

This postulated characterization is attained through the construction of two geomet- 
rically meaningful operators on the phase space C of these equations: an almost-Poisson 
operator /C, and a recursion operator 1Z which can be both realized in terms of C+ and 
C- (see l|44|) and jSJ ), and as K,J~ X (where J — Diag(— is the standard Poisson 
operator on C.) 

The former depiction of 1Z 7 allows us to view the 7?.-generated hierarchy (|65|l as a sub- 
hierarchy of the Chiu-Ladik flows. The latter realization of this recursion operator, attests 
to its true geometric character, and it is responsible for elucidating the bi-Hamiltonian 
character of the hierarchy. 

The skew operator K. can be portrayed as the discrete counterpart of Magri's Poisson 
operator for NLS, based both on its mode of construction and on its geometric properties. 
For example, we will see that, once their bi-Hamiltonian nature is established, the skewness 
of /C serves to prove the commutativity of the AL flows generated by 1Z. 

4.1 The Geometric Context: AL as a Hamiltonian System 

We begin by defining the phase space C of the AL equations as the space of complex, 
vector valued functions of a discrete real variable k G Z given by 

C = {{q k ,r k ) T e ••• x C 2 x C 2 x C 2 x ••• : \q k \,\rk\ -> as k ~> ±oo, }. 

A point in C will either be denoted by (<?, r) T = [■ ■ ■ , (q k ,r k ) T , (q k+1 ,r k+1 ) T , • • • ] , or just 
by (qk,i"k) T , depending on the context. 
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As in the continuous case, we may endow the tangent-bundle to C, with a non- 
degenerate bilinear form (•,•) arising from a point-dependent inner product defined by 

((a, bf, (c, d) T ) (qr)T = - ^ (a k d k + b k c k ) (1 - q k r k )-\ (60) 

on the tangent space to C at (q,r) T ; that is, (a,b) T , (c,d) T lie in T^ ^tC ~ C. Under the 
identification r k = — q k , (|60fl defines a positive-definite, real inner product. 

Given a functional F over C we may define the (discrete) functional gradient of F, V-F, 
hy VF((q,r) T ) = [• • • , W k F((q, r) T ), V k+1 F((q, r) T ), ■•■], where for k in Z, 

V fc F((g,r) T ) = -(l-q k r k )a 1 (5F/5q k ,6F/5r k ) T = VF(( % , r fe ) T ), (61) 

where ci is the usual Pauli matrix with ones in the off-diagonal. 

Now, let Jtq r)T = Diag(- • ■ , J k ,Jk+i, ■ ■ • ), Jk = Diag(— i, i) for all k, define the (stan- 
dard, point-independent) skew-symmetric operator on Tr q>r \TC, acting as multiplication by 
J k on each C 2 component. The operator J (which may also be regarded as an operator 
on the tangent bundle to C) may be used to define the following Poisson bracket on C: 

{F,G}j((q,r) T ) = (VF,JVG)((g,r) T ) 

= ((l-%r fc )(^/<5r fc ,(5^/%) T ,Diag(-i,i)(l- ( 7 fc r fe )( ( 5F/ ( 5r fe , ( 5F/ ( 5 (Zfc ) T ) (9i , )T , (62) 

where, F and G are smooth functionals over C which become real- valued under the identi- 
fication r k = — q k . Note that as in the continuous setting, this reality requirement implies 
that F and G must be either symmetric or anti-symmetric in r k , q k (for all k.) 

Remark 4 Since the operator J k defined above does not actually depend on k, we will 
often use the notation J to refer to the 2x2 complex matrix Diag(—i,i) originally de- 
fined as J k - By the same token, from now own we may omit noting explicitly the point 
dependence of operators which are clearly characterized as such by their definitions. 

Let H((qk,r k ) T ) — — X^oofe^+i + r k q k+ i). Then the AL equations (|36|l can be ob- 
tained as a Hamiltonian system on C with respect the ^7-bracket (I62|l for the Hamiltonian 
functional H. In this same manner, the Hamiltonian H((q k ,r k ) T ) = —^^^iqkrk+i + 
rkqk+i + 21og(l — qkrk)) produces the standard (vector) AL flow, which becomes in 
the reduction r k = —q k . (See also [TT].) 
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4.2 The Operators K, and 7Z and the AL Hierarchy 



Guided by the parallel with our continuous setup and the form of the operators C + and 
we re-consider the operator 1Z defined in (|52[l . and define a new operator 



fC((Qk , r k ) T ) = TZ((q, r) T )J = -i{ 2(E+ + ET) 



+ (1 - r k q k )(E+ - ET) 



1k{J k - J k ) (l^-) -»( J l + - J WIT^ ( 
M4 Jk){T^) -rkiJt-J^l^ 



/ qk r k ~1k 

+ (l-r k q k ){E + + E-)[ 1 -lf* 



9k(Jjt-Jk)(rj) q k (J+ -J k ){qj) 
-r k {4 - J k ) fa) -n(J+ - J fc ~) 



(£+-£T) 

(63) 



Just as with 7Z, K. is point-dependent and non-local; its definition at site k involves func- 
tions of the vector potential at neighboring sites, for example, functions of q k ±i, r k ±i- 
As /C acts on the tangent bundle of C, it may be used to define the bracket 

{F,GM(<z fe ,r fe ) T ) = <VF,/CVG)((<?,r) T ), (64) 

on the class of functionals over C just discussed. 

Theorem 3 The operator K, is skew-symmetric. 

Proof. One first shows that 1Z* — o 3 lZa 3 where, as before, * denotes the adjoint of an 
operator with respect to the inner product jBDjl. and 03 = Diag(l, —1). This calculation 
uses the fact that both ±102) and 03 are self-adjoint, while C~ = a^C^a^ 1 . 

The result now follows from the definition of K, in terms of 1Z and J = —io~%. Indeed, 

{F,G} K = (VF./CVG) = (VF,1ZJVG) = {1Z*VF,JVG) = (-J(a 3 TZa 3 )VF,VG) 

= (lZ(icr 3 )VF,\7G) = -(TZJVF,VG) = -(VG,/CVF) = -{G,F} K . □ 



Corollary 2 The IC-bracket (|64f) defines an almost Poisson structure on C . 

Proof. In addition to being skew, the K, — bracket satisfies the Leibnitz identity, as we 
show in 22]. □ 
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4.3 The bi-Hamiltonian character of the AL equations 

Let us now consider the functional H ((q k , r k) T ) = — E^oo l°s(l ~ Qk r k), on C, and 
observe that JVH ((q kl r k ) T ) = i(q k , —r k ) T '. If we now define Xr n \ = lZ n J\7H , then 
Xi n \ defines a hierarchy of fields on the tangent bundle to C which, we claim, has the 
following properties: 

(i) The fields Xi n ) mutually commute; 

(ii) The hierarchy comprises the AL equations l|36|) . which occur as Xny, 
(Hi) The evolution equations 

((q k ,r k ) T ) t = K n JV k H = X (n) , (65) 

are equivalent to a sub-hierarchy of the family l|45(l and are therefore integrable; 

(iv) The fields X^ in the hierarchy are Hamiltonian with respect to J and also with 
respect to K,, i.e., bi-Hamiltonian. 

Remark 5 As indicated in the term bi-Hamiltonian here refers to sequences of fields 
X( n ) which satisfies the Lenard relations X/ n -\ = /CVi/„_i = JVH n '91. The clarification 
is pertinent due to the fact that the Poisson nature of the skew operator K. has not yet been 
established. 

It is a good exercise to verify items (Hi) and (iv) explicitly for the first 2 iterations 
[T%| (in the process, item (ii) becomes apparent.) As such calculations indicate, being able 
to show that the fields £™ ((r k , q k ) T ) are in the kernel of C+ — C- for all n is of central 
importance in proving claim (Hi). This is precisely the content of Theorem 0] In it, we 
use once again the identities derived in Appendix II to extend the Calogero-Degasperis 
"kernel-condition" result (see ©) to the discrete setting. CorollarylSl the discrete-analog 
of Proposition [3 will establish item (Hi) in full generality. Theorem will then use the 
resolvent identities l|53[) . I|54|) together with Corollary [3] to establish (iv). The proof of 
item (i) is rather immediate given (iv). It is the content of Corollary [5] 

Throughout the remainder of this note, E /,• = Efc-00' a ^ so ^ ne difference operator 
£+ — C- will often be written as £+ — C- =T>\ +D2, where 

Z>i-M r,m)(F+ F-\f r kT,Qj/(^-rjqj) -r k E »v/(l - r 3 q 3 )\ 

T>\ - (1 r k q k )(E * ) (_ ftE<&/( i _ P ^) qk ^ rj /(l - r j9j ) ) ' (66) 

and 

V 2 = ( Tk p qj Tk p r ^ (E+ - E-). (67) 

Lemma 1 Let j = 1,2. IfVj (L m (r k ,q k ) T ) = Vj ((L- 1 )" 1 (r k , q k ) T ) = for all m, then 
(£+ — £_) (£™ (rfc, qk) T ) = for all n (m, n non-negative integers.) 
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Proof. This follows immediately from the definition of £+ (see 11441 above,) the inverse 
relation between L and L^ 1 (described in |18| . proven in |S]) and the binomial theorem. 
□ 

Proposition 5 For V l as in V l (L m {r k ,q k ) T ) = V 1 ({I,- 1 )" 1 {r kl q k ) T ) = 0. 

Proof. Recall that the vector of squared eigenfunctions ((3 k , —a k ) T defined in Appendix 
II may be explicitly viewed a generating function for the hierarchy of fields arising from 
powers of L (through the resolvent identities in Theorem[2]) The same is true of (—j3 k , a k ) T 
and the fields given by L _1 ((rfe, q k ) T ). Using the inner product ((60(1 one observes certain 
relations between the aforementioned fields, which can be compactly expressed as: 

((1 - r k q k ){p k , -a k ) T , -ia 3 (l - r k q k )(f3 k , -a k ) T ) = 1)^(1 - r k q k )((3 k a k - a k (3 k ) = 0, 

^ (68) 
((1 - r k q k )(-pk, a k f, - r k q k )(-p k ,a k ) T ) = 0. (69) 

We note that the previous identities translate into commutativity relations among special 
sub-classes of Hamiltonian functions, namely {H^, H^} = {H!^,H®} — (refer to ^B] 
for the precise definitions of and H®.) In addition, one may focus attention on the 
particular relations 

= ((r k ,q k ) T , -ia 3 (l - r k q k )(j3 k , -a k ) T ) = i^(q k f3 k + r k a k ), (70) 
= ((r k ,q k ) T , -10-3(1 - r k q k )(-/3 k ,a k ) T ) = -i^(q k (3 k + r k a k ), (71) 

obtained from expanding the first entry on the left-hand-side of ((68(1 (respectively ((69(1 ) 
in powers of z~ 2 (respectively z 2 ) as suggested by the resolvent formula (|53|l (respectively 

Equation 1)70(1 and ((66(1 together with ((71(1 and 1(66(1 immediately show that 

V x ((1 - r k q k )(p k , -a k ) T ) = 0, and V 1 ((1 - r k q k )(-0 k , a fc ) T ) = 0. (72) 

Using now ((53(1 and the first part of ((72(1 . we see that 

= V 1 (r k , q k ) T + z- 2 V x (P(r k , q k ) T ) +■■■ + z^V, (L m (r k ,q k ) T ) +■■■ 

or, equivalent^, V x (L m (r k ,q k ) T ) = 0, for all m. The identity V x {{L- x y n {r k ,q k ) T ) = 
stems from Eqs. 1(54(1 and the second piece of 1(72(1 in a parallel manner. □ 

Theorem 4 For any non-negative integer n, (£ + — C-) (£™ (rfc, gx) T ) = 0. 
Proof. Due to Lemma ^ and Proposition [5] it suffices to prove that 

V 2 ((L- 1 ) m (r k ,q k )) = V 2 (L m {r kl q k )) ( => 0. (73) 
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Below, we present the argument showing that equality (b) in l(75|) holds. An analogous 
argument (based on the identities given at the end of Appendix II) shows that identity 
(a) holds as well, validating the theorem. 

The key step in the argument closely resembles the recurrence strategy employed in the 
derivation of the resolvent identities 153|) . (|54|l . as described in Appendix II. In the current 
scenario, the technique entails re- writing the vector (E + — E~) [(1 — r k q k ){fi kl — «fc) T ] as 

(z- 2 - z 2 )Diag(l, -1) [(1 - r k q k )(l3 k , -a k f - [r k ,q k ) T ] (74) 
+Diag(r fc+1 + z~ 2 r k , q k+1 + z 2 q k ) (q k f3 k + zj k + 1, -(r k a k + z~ 1 5 k - 1)) T (75) 
-Diag(-(r fe _i + z 2 r k ) 1 q k _ 1 + z~ 2 q k ) E~ (q k (3 k + zj k , r k a k + z~ 1 6 k ) T , (76) 

by considering the difference between equations H1(J(J[) and IjlOlU in Appendix II. 

The remainder of the proof amounts to a series of calculations. First we utilize Eqs. 
(|74[1 - (|76[) together with (|98|l and 1)1 U4|) in Appendix II, to obtain a simplified three-term 
expression for X>2((1 — r k q k )((3 k , —a k ) T ). We then exploit the particular form of the re- 
sulting expressions to show that their combined sum vanishes. Finally, the same argument 
used at the end of the proof of Proposition shows that the vanishing of the vector 
2?2((1 - r k q k )((3 kl -a k ) T ) implies the vanishing of T> 2 (L m (r k ,q k ) T ) for all m. 

Due to (|70|l . the matrix pre- factor of the operator T>2 applied to expression l|74|l . 
becomes: 



-r k q k 



-a k 



(77) 



o qtc) \£qj -Er?, 

Expression l|77|) is the first term in the sought-after three-term expression for T>2((1 ~ 
r k q k )(fi kl —a k ) T ). The second and third terms (appearing in J7SJ, (dJ) are essentially 
obtained by re- writing l|75[) and i|76|) by means of i|104l) . in Appendix II. Specifically, the 
matrix pre- factor of T>2 applied to expression (|75|l . yields the new version of the second 
term in the expression for £>2((1 — r k q k )(f3 k , —a k ) T ), namely, 



'7,- 



-l+z 2 r k 









q k+ i + z 2 q k 



■z-\){l,l 



-a k 



qk) VE* J2 r j 
- J fc+lfe) 

=Diag(r fc , q k ) qj(r j+ i 

+ [~1j(r J+ i + z~ 2 rj) + rj{q ]+1 + z 2 ^-)] J+ +1 {qi(3i + naj)(l, 1) T 
=Diag(r fe , q k ) ^ [q 3 {r 3+ i + z~ 2 rj) 

+ {-Qj-i{ r j + z ~ 2 rj-i) + rj-ife + 2 2 %-i)) J/(*A + nan)] (1, if 
=Diag(r fe , q k ) ^ [q 3 {r 3+ i + z~\j) 
+ (ijPj + r 3 a j)J7 (-qi-m + n-iqi 



(78) 
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where, recall, J^(uj) = Ylj=-oo u j- The last equality in the above string of identities 
results from a direct application of the summation by parts formula (see Appendix A in 
|18|.) The previous one, amounts simply to a convenient shift in the summation index for 
the bi-infinite sum. This same series of calculations leads to the desired third and final 
term. That is, using I|1U4|) . and summation by parts, the action of the matrix pre- factor 
of T>2 on l|7bj) leads to the vector 

Diag(r fc , q k ) ^2 [ ~ r j(Qj-i + z ~ 2( lj) 

+ (fcflj + rjOj)Jj (-n-i« + qi-m + (z- 2 - z 2 )qin) ] (1, l) T (79) 

Consider now the sum of expressions (|78() and (|79() . The J~ -independent terms all 
disappear: the two z-dependent ones cancel out; the remaining two telescope away, leav- 
ing no boundary terms. Focusing now on the ./"-dependent terms, one observes that 
the ^-independent ones cancel off identically. The four z-dependent terms are pairwise- 
telescoping, and hence reduce to Diag(rfe, q k ) [(lj0j+ r j a j)lj r j( z ~ 2 ~ z2 )] (1j 1) T ' which 
is exactly the negative of i|77|) . ft follows that, for all (non-negative) m, 

V 2 ((1 - r k q k )(p k , ~a k ) T ) = 0, and hence, V 2 (i m (r fe , q k f) = 0. (80) 

A parallel argument, based on identities (|106J) . (|107|1 . (|108|1 in Appendix 11 yields 

V 2 ((l~r k q k )(~$ k ,a k ) T ) = 0, and hence, V 2 (L- im (r k , q k ) T ) = 0. (81) 

The theorem follows. □ 

Corollary 3 // (£+ - £_) (C^{r k ,q K ) T ) = for m > 0, then 

(_°! ~o) JWkH ° =2n (-l J)^^ ~^JV k H , (82) 

for n>0 and H = - Sfcl-oo ^Si 1 ~ r fc9fc)- 

Proof. The statement may be verified by direct computation for n = 1. One may then 
proceed by induction on n. Assuming that identity i|82[l holds for n — 1, we show below 
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that it also holds for n. Theorem 0] is used in the last line below. Indeed, 



1 

-1 



1 

-1 



(£+ + £-) 
(£+ + £ 



1 



JV k H = i 



1 

-1 



-1\ / 1 



i (i ; v-i o ' + 

„_! / 1\ n _i / r k 



(£++£_)(£++£_)"■ 

i f r*. 



-1 / rfc 



9ft 



i2 



-1 



i2 r ' 



iT 



1 

-1 

1 

-1 



= i2 r ' 



1 
-1 

1 

-1 



CI 



Qk 

-1 

1 



-n-1 I r k 

qk 



□ 



Corollary 4 (Validity of Item (Hi).) The evolution equations (|65|l are equivalent to a 
sub-hierarchy of the Chiu-Ladik family (|45|l and, thus, integrable. 



Proof. Working from Eq. (|65|l and using Corollary|21 we have: 



1 



K n JV k H Q = il_ 1 o ){C++C-T ( I 



-1 



qk 

-Tk 



1 



1 +\q k 



or, equivalently, ((~r k ,q k ) T ) t — *2" 'Cl(r k ,q k ) 7 ', which is of form l|45|l . as claimed. □ 

Theorem 5 (Validity of Item (iv).) The fields -X7 n ) in the hierarchy <|65|) are Hamiltonian 
with respect to both, J and K,, i.e., bi- Hamiltonian in the sense of Remark\^ 

Proof. Let us once again consider the resolvent formulas for L and Starting from 

(1531) and matching inverse powers of z 2 , we see that: 



(1 - r k q k )6 k H Q = (r k ,q k ) T , and also, (1 - r k q k )6 k H" = L j (r k , qk) 1 
Using the definition of given in (ltj 1 1) - we observe 



(83) 



JV k m =Diag(-M) 



1 

-1 



I? 



1 

1 

-1 

1 



1\ q k 

1 \r k 



(Ik 

-Tk 



1 

-1 



V 



-1 

1 



Reasoning in a parallel fashion for the L 1 resolvent, (|54|l . we get: 



(1 - r k q k )5 k H] = (L~ ) 



o 1 

-i o 



-1 

1 



JV k H Q . (84) 



JV k H . (85) 



2G 



(Notice that, based on Eqs. and J57JI, one has S k Hf = Sf-Cj, S k Hf = 8 k C 3 , for j > 0.) 
Based in Eqs. (|8"^|) and (JHSJ, we have: 



= TZ n JW k H c = E1 2 r 



= 2 r ' 



'fc-no 

1 

-1 

1 

-1 



1 

-1 



-1 

1 



JV k H 



E 



-1 

1 



i — k>k 



2 n JV k 



i — k>k 



jja 

n n-2k 



JV k H 

L ( n -k)-k + 

E 



n—k<k 



n — k<k 



jja 

n -(n-2k) 



-1 \ k—{n — k) 



JV k H n 



-1 

1 



and so, the last equality, explicitly depicts the alluded hierarchy of flows as j7-Hamiltonian. 
Given this, one uses the fact that 1Z = K,J _1 in a recursive fashion to show that these 
flows are also /C-Hamiltonian (in the sense of Rcmark[SJ) Indeed, for n = 1 and n = 2, 

njV k H = JCV k H = JV k H u and -R 2 JV k H a =JCJ- 1 K,V k H Q = lCV k H 1 , 

respectively. It then follows via an inductive argument, that lZ n J"V k Ho = K,V k H n -i = 
JV k H n , for n > 0, where H = H . □ 



Corollary 5 (Validity of Item (%).) Let J and K. be skew operators on the tangent bundle 
to C, and let {Hj} be a sequence of Junctionals on C indexed by j in Z. If J^V H k+ i = 
K,VH k , then the Hj pairwise commute, that is {Hj, H k }j = 0. 

Proof. Simply using the skew-symmetry of K. and J , we have: 

{H„H k }j = (VH 3 ,JVH k ) = (V^./CViffe.x) - - (KS/Hj, ViT fe _i) 

= - {JVH j+ i,VH k -i) = (VH J+1 ,JVH k ^) = {H j+1 ,H k _ x }j. (86) 

Assume now k > j. Then after {k — j) iterations of the procedure yielding (|86[l . we obtain 
{Hj, H k } j = {H k , Hj}j, which implies {Hj, H k } j — 0, as the ./-bracket is skew. □ 
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Appendix I: Key elements in the proof of Theorem [T] 

Let (<ik,rk) T = ((Xfe, bk) T denote an arbitrary variation of rfe), and ^fc denote a (vector) 
solution to the eigenvalue problem (|18(l . The induced variation on i/^ satisfies 

(87) 

Setting z>fe = $(fc)tife for $(fc) as in (|23H . and using the standard variation of constants 
method to solve (|87|) for yields: 



Upon inverting — 1), one obtains 



u „ = - £ $(fc + l)- 1 (? ^W + c, !)„ = $(«)«„. (89) 

k—n ' 

In particular, choosing v n = ?/;„, we see that as n — > oo, "0 n ~ (0, z~") T . So ■0 n vanishes 
and, clearly, so does the sum in l|89l) . which implies c = 0. Multiplying the resulting 
equation through by z n and focusing on the variation of just the second component of the 

(2) 

resulting vector (i.e., z n ipn ') we see that: 



1(2) 

(90) 



^[(z"$ (n)) $ (fc + i)-i^ 

OO 



where 



(91) 

,J k+i ^k+i J 

(Here, one uses relations Q to determine det $(fc) = (j^^f ] -tfV^ = a/P(k) [TBI EI) 
Taking the limit of (|90|l as n — ► — oo and using the asymptotics (|23|) together with (|91l) 
above, we get: 

oo 

(a/Co) = £ (P(fc + 1)/C ) r k + 4+>£ 2) **) ■ (92) 

k— — oo 

Multiplying by Co/a, we obtain (@U and, hence, (gTJ. 
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Focusing on the variation of z n $n and arguing in a similar manner, one arrives at 
the identity 

oo 

(a/C )= ( p ( k + !)/ C o) (^IA 1] h ~ ffii$ 2) 9k) , (93) 

k— — oo 

Theorem ^ follows. 

Appendix II: Key elements in the proofs of Theorem [2] 
and Theorem |U 

We begin by defining the following matrices of squared eigenfunctions: 

r, -< T „ i° Mr'-K'f -<P ( M ) \m_(c k -V, 
v "-*"\p -i)*k -[^Ur -JtWJ a V s * -D*)' (9 ) 

Wk -^{o -i)*^-[^¥ k i -4 2) 4ij a -u (95) 

Note that 14 's entries are standard squared eigenfunctions, whereas W k is composed of 
semi-shifted products. Observe also that the constant matrix defining these identities is 
chosen so that the off-diagonal entries of W k are precisely the components of the vector of 
squared eigenfunctions defining the generating gradient S k log a + 6 k Ho. 
One can directly verify the following identities for V k and W k ' 

• W k+1 = £ k W k £^ 1 ; (96) 

• W k = V&k*^ = V k S?; (97) 

• (r k a k + Z - l 5 k ) - {q k fJ k + z lk ) = 1, (98) 

where £ k is, recall, the coefficient matrix for our eigenvalue problem (|18|) . and the entries 
of W k are now written in terms of a k , f3 k , jk and S k , as in (|95|l . 

By writing (|96|l in terms of the new form of W k and tracking only the off-diagonal 
entries of the resulting matrix, one observes that 



(1 - r k q k ) 



~OL k 



r 2 0\ / /3 fc \ fr k r k+1 a k + z 1 (r k j k + r k+1 S k ) 
z 2 J \-atk I I -9k1k+xPk ~ z(q k S k + (?fe + i7fc) 



(99) 

Adding and subtracting the expression Diag(z 2 , z 2 )[—r k q k (/3 k , —a k ) T — (r k , q k ) T ] on the 



29 



right-hand-side of (|99|l . and using identity Q98[l. one may re- write it as 

- 2 



(1 - r k q k ) 



Pk 
-oi k 







z 



(1 - ^feft) 



Pk 

~Ot k 



(100) 



rfe+i + z l r k 







ft/?fc + Z7 fc + 1 
-{rkctk + z~ 1 6 k - 1) 



By multiplying i|100|) by Diag(z 2 , z 
tion in terms of the opposite shift: 



q k+ i + z z q k/ 

~ 2 )E~ , we obtain an almost perfectly symmetric equa- 



E' 



(1 - r k q k ) 



Pk 
-a k 



z 2 
z~ 2 

r k -i + z' z r k 




(1 - r k q k ) 




Pk 
-oik 



(101) 



q k -i + z 2 q k 



E' 



-(qkPk + zik) 
r k a k + z^Sk 



Spelling out the relationship given in l|97|l one discovers: 

D k = z~ 1 5 k - qkPk, and - C k = r k a k - zj k , 



(102) 



so that (|9"8|) may be written as: D k — C k = (z~ 1 5 k — q k fi k ) — {r k a k — zj k ) = 1- Using 
this fact together with the previous relations (|102|) . we see that the dry-dependent factors 
of IjlOOl) and l|101|) may be written as 



(qkPk + r k a k + C k ) 



1 

-1 



and (q k Pk + r k a k + C k ) L + 



(103) 



respectively. Finally, writing out V k +i in terms of its entries, we find that C k satisfies a 
difference equation given in terms of a k and f3 k , namely 

{E+ -l)C k = q k p k +r k a k . 

Taking C k = —^JLkljPj + r j a ji and substituting this expression into Eqs. (|103f) . we 
obtain 



qkPk + z'fk + 1 
-(r k a k + z^5 k - 1) 



q k Pk + zjk \ _ f-J k+ i{<lj) J k+ i( r j) 



r k a k + z 1 4 



~ J k+M) J k+i( r i)J \- a k 



Pk 



Pk 
-oi k 





(104) 



expressing the ^7-dependent factors in l|100fl and (|101() in terms of the sum operator j£ 
in the formula for L, and the non-standard squared eigenfunctions a k , [3 k . 

Equations l|100l) . I|101|l used in conjunction with (|104|l (as well as (|10uf) . (|107|l together 
with l|108|) ') play key roles in the proofs of Theorem|2|and Theorem^] as indicated therein. 
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The latter triple arises by considering the matrix of (semi-shifted) squared eigenfunctions 
tied to the generating function (—/3 k ,a k ) T described in (|43|l . namely, 



1 - ?(2) :f 2 y 7(2) W , (n^fc + z 4) - (g fc /3 fc + z lk ) = 1, 



-0, 



(105) 

Starting with (|105|) . and applying a procedure analogous to the one described in this 
appendix, we obtain: 



(1 - r k q k ) 



-fa 
&k 



z~ A 
z 2 



(1 - r k q k ) 



-fa 
a k 



r k +i + z 2 r k 







q k+1 + z*q k 



qA + z% + 1 
'{r k a k + z _1 5fc - 1) 



(106) 



(1 - ntft) 



-ft 



z z 

Z- 2 



(1 - r k q k ) 



-fa 
a k 



r k _i + z A r k 




gjfc-l + z 2 q k 



qkfa + z% 
~(r k a k + z^ x l k )) ' 



where 

f qkfa + z% 
V-(^afc + ^ _1 4) + l 
as just mentioned above. 



-Jk+i(rj)\ (-fa 



(107) 
(108) 



5 Conclusion 

To summarize, this paper shows that the AL hierarchy can be explicitly viewed as a hier- 
archy of commuting flows which: (a) are Hamiltonian with respect to both the standard, 
local Poisson operator J, and a new non-local, skew, almost Poisson operator JC, on the ap- 
propriate space; (b) can be recursively generated from the recursion operator 1Z = K,J~ X . 
In addition, the proof of these facts relies upon two new pivotal resolvent identities which 
suggest a general method for uncovering bi-Hamiltonian structures for other families of 
discrete, integrable equations. 

Another result stemming from the current research is the clarification of the geometric 
framework that underlies a certain class of geodesic linkages evolving on the sphere |10U18j . 
A linkage on a Riemannian manifold is essentially defined by specifying a sequence of points 
connected by geodesic arcs. A closed linkage is usually called a polygon. Such linkages 
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are related to the AL hierarchy via the evolution for their "discrete" geodesic curvature 
10 . In this regard, Lozano's preliminary results include a geometric interpretation of 
a compatibility condition associated to a Lax pair for AL and, consequently, a bijective 
correspondence between discrete, integrable mKdV flows (also AL flows) and linkage flows. 
(For details on this, see a ls° see |lUH15llT7j for background in terms of continuous- 
analogs of the linkage models.) 

Let us now close by summarizing some of the many possible avenues for further research. 
First off, a definite answer to the question of whether or not /C actually defines a Poisson 
structure would be desirable on several counts. If K. were Poisson, one could turn to 
exploring the possible connections between the /C-induced Poisson bracket and the bi- 
Hamiltonian structure for finite AL described by Faybusovich and Gekhtman ^21- This 
work considers the AL hierarchy within the larger class of full Toda flows in sl(n) and 
presents a bi-Hamiltonian formulation for (finite) AL stemming from the bi-Hamiltonian 
structure of these Toda flows. 

On the other hand, an obstruction to the Jacobi identity would place the /C-bracket in 
the category of almost Poisson structures and could perhaps steer the investigation in the 
direction of non-holonomic mechanical systems [5]. Such systems possess an underlying 
Hamiltonian structure that is (strictly) almost Poisson. 

The elucidation of new bi-Hamiltonian structures and their connection with the evo- 
lution of non-stretching classes of linkages could also be pursued in the context of the de- 
focusing Ablowitz-Ladik system (obtained from in the reduction r k =q k ) and other 
discrete integrable equations. This would serve as a test of the robutness of our methods 
for deriving operators such as K,, 1Z, and resolvent identities such as those obtained for 
L and L~ l . It would be good, for instance, to understand how recursion operators (such 
as 1Z, or even L and L~~ r ) are encoded in the squared eigenfunctions of a linear problem 
(through resolvent identities of the right kind.) 

The connection with "physical linkage" spaces could also be pursued further, both 
for discrete, integrable mKdV, AL and potentially for other discrete integrable equations. 
In the context of discrete integrable mKdV, AL and non-stretching spherical linkages 
for instance, one should aim at understanding the linkage recursion schemes proposed 
in JSj m Poisson-geometric terms. Ideally, a well-defined lift of TZ to the space of non- 
stretching linkages could be defined and then parsed out as the composition of two Poisson 
(or perhaps one Poisson and one almost Poisson) operators. Connections with results of 
Langer and Pcrlinc (for the case of continuous NLS and the FM model US])) and Kapovich 
and Millson (regarding the symplectic geometry of non-stretching polygons |13H14j . could 
also be explored and addressed in the appropriate context. 
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